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Abstract 

We study the decay of Df mesons into final states involving a Kg and report the 
discovery of Cabibbo suppressed decay modes Df — ► Kgir~ tt + tt + (179±36 events) 
and Df — > Kgir + (113±26 events). The branching ratios for the new modes are 

rg^g;-;^) = 0-18±0.04±0.05 and j^ggff^ = 0.104±0.024±0.013. 

PACS numbers: 13.25.Ft, 14.40Lb 



An essential ingredient to accurately model backgrounds in heavy quark sys- 
tems involves the identification and categorization of missing decay channels 
in the charm sector. This is particularly important for the Df decays where a 
substantial part of its hadronic decay rate is yet to be identified. Only two Df 
Cabibbo suppressed decays have been reported, namely Df — > K + 7i + n^ [1,2] 
and its resonance substructure and Df — > K + K + K~ [3]. It was found that 
°t^*tl +n ~ + = 0.127 ± 0.007 ± 0.014 and %^ K ^ K \ = (8.95 ± 2.12 
— liii) x 10~ 3 . The two Cabibbo suppressed channels differ by an order of magni- 
tude (partly due to phase space) and additional decays are needed to establish 
patterns. In this paper we report the discovery of two Cabibbo suppressed de- 
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cays of the meson; — > K s n 7r + 7r + and Df — > K s n + . No inclusive esti- 
mates of the branching fraction for Df — > K s tt~7t + tt + have been reported, but 

several predictions exist for the branching ratio of ^+^^0^+ [4,5,6]. Through- 
out this paper, charge conjugate modes are implied unless explicitly stated 
otherwise. 

II. THE FOCUS EXPERIMENT 

The data come from 6 billion events recorded during the 1996-1997 fixed tar- 
get run at Fermilab. Electrons and positrons with an endpoint energy of ap- 
proximately 300 GeV bremsstrahlung, yielding photons which interact in a 
segmented beryllium-oxide target to produce charmed particles. The average 
photon energy for events which satisfy our trigger is approximately 175 GeV. 
Charged particles are tracked and momentum analyzed by a system of silicon 
vertex detectors [7] in the target region, multi-wire proportional chambers 
downstream of the interaction region, and two oppositely polarized dipole 
magnets. Particle identification is performed by three threshold Cerenkov 
counters, two electromagnetic calorimeters, a hadronic calorimeter, and two 
muon systems. The main FOCUS trigger required tracks outside of the central 
region and approximately 25 GeV (or more) of energy in the hadron calorime- 
ter. 

Df decays are reconstructed using a candidate driven vertex algorithm [8]. 
A decay vertex is formed from the reconstructed charged tracks. The K s — ■> 
7r + 7r~ decays are reconstructed using techniques described elsewhere [9]. Briefly, 
K s — > 7r + 7r~ decays can occur anywhere along the spectrometer. Depending 
on where the decays occur (upstream of the first magnet or inside the magnetic 
field) and on how many multi-wire proportional chambers each pion passes, 
the Kg are given a type number and the different types vary in mass resolu- 
tion and in purity. The momentum information from the K s and the charged 
tracks is used to form a candidate D momentum vector, which is intersected 
with other tracks to find the primary (production) vertex. Even though it is 
possible for the production vertex to be identified with a single track plus 
the Df momentum vector, the signal quality is greatly improved by demand- 
ing at least two primary tracks. Events are selected based on several criteria. 
The confidence level for the production vertex and for the charm decay ver- 
tex must be greater than 1%. The likelihood for each charged particle to be 
a proton, kaon, pion, or electron based on Cerenkov particle identification is 
used to make additional requirements [10]. We define a x 2 -like variable Wi as 
— 2 ln(likelihoodj) for the hypothesis i. In order to reduce background due to 
secondary interactions of particles from the production vertex, we require the 
decay vertex to be located outside the target material. We enhance the sig- 
nal quality by cutting on the isolation variables, Isol and Iso2. The isolation 
variable Isol requires that the tracks forming the D candidate vertex have a 
confidence level smaller than the cut to form a vertex with the tracks from 
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the primary vertex. The isolation variable Iso2 requires that the tracks not 
assigned to the primary or secondary vertices have a confidence level smaller 
than the cut to form a vertex with the D candidate daughters. 

III. D+ -> K° s 7i-7i + 7i + CHANNEL 

For this channel we have excellent secondary vertex resolution with at least 
three charged tracks defining the vertex. We require Iso2 less than 1% so 
the secondary vertex is isolated from other tracks. We require I sol less than 
1% to make sure the Df tracks do not originate at the primary vertex. The 
reconstructed mass of the Kg must be within four standard deviations of the 
nominal K s mass. The typical Kg mass resolution is approximately 6 MeV/c 2 . 
For each pion candidate we require a loose cut that no alternative hypothesis 
is greatly favored over the pion hypothesis: min (W e , W K) W p ) — > —5. 
For the charged kaon candidate in the normalization channel we require W w — 
Wk > 2. We also require the distance L (~5 mm) between the primary and 
secondary vertices divided by its error (~500 jum) to be at least 7. Lastly we 
require an additional D* + —D° cut for the K s 7i'n + n + sample. The K s n^n + n + 
invariant mass minus the highest K s ii~ii + mass combination must be greater 
than 0.160 GeV/c 2 . This eliminates D* + background events, which simplifies 
the fitting function. 

Figure 1(a) presents the invariant mass plot for the normalization channel 
K s K^n + n + which is the cleanest four body decay containing a Kg. The 
figure contains the Cabibbo suppressed channel from the D + as well as the 
Cabibbo favored signal. We fit the D + and Df signals with Gaussians. 
We include a background contribution from D + — > K s 7i^ii + 7i + where the 7r~ 
is misidentified as a kaon and the shape is determined from a Monte Carlo 
simulation. The combinatoric background is fit with a 2 nd degree polynomial. 
We find 763 ± 32 Df signal events at L/a L > 7. It is worth noting that 
this channel has been previously studied by the FOCUS Collaboration and 
the signal yields reported in this paper are comparable to the results already 
published [11]. 

Figure 1(b) shows the K s n~-K + iT + invariant mass plot for events that satisfy 
the above cuts. The plot is dominated by the Cabibbo favored decay D + — > 
K s it~7t + 7t + while the signal is barely visible. Figure 1(c) is the same 
Kg7r~7T + 7T + invariant mass distribution in the region above the D + peak. The 
Figure 1(c) mass distribution is fit with a Gaussian with the width fixed from 
Monte Carlo for the signal and a first degree polynomial for the background. 
A signal of 179 ± 36 D+ events is found from the fit. 

We measure the branching fraction of the — > Kgii^7i + 7i + mode relative 
to — > KgK~7T + 'n- + . The relative efficiency is determined by Monte Carlo 
simulation. The relative branching fraction is reported assuming non-resonant 
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Fig. 1. Invariant mass distributions for (a) KgK~ir + ir + (background reflection from 
a mismeasured pion from D + — » Kgir~ir + 7r + is included), both the D + and -D+ 
signals are evident, (b) K^ir~TT + ir + (not fitted to show the large Cabibbo favored 
D + contribution), (c) K^tt~-k + it + (with the invariant mass only plotted above the 
D + mass). The mass distribution is fit with a Gaussian with the width fixed from 
Monte Carlo for the signal and a first degree polynomial for the background. 



decays for both channels. We test for dependency on cut selection in both 
modes by individually varying each cut. In Figure 2 we present the ratio of 
branching fractions for — > Kg7i^7i + ii + relative to — > KgK~7r + ir + as 
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a function of significance of separation between the primary and secondary, 
isolation of the secondary, and confidence level of the secondary vertex. 
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2. The ratio of branching fractions for D s 



Kqir 7r + 7r + relative to 



Fig. 

Df — ► XgA' _ 7r + 7r + as a function of significance of separation between the pri- 
mary and secondary (first eight sets), isolation of the secondary (next five sets), 
and confidence level of the secondary vertex (final 10 sets). 



We studied systematic effects due to uncertainties in the reconstruction effi- 
ciency, in the unknown resonant substructure, and in the fitting procedure. To 
determine the systematic error due to the reconstruction efficiency we follow a 
procedure based on the S-factor method used by the Particle Data Group [12]. 
For each mode we split the data sample into two independent subsamples 
based on momentum, particle versus antiparticle, decays inside the tar- 
get material versus outside of target material, and on the period of time in 
which the data was collected. These splits provide a check on the Monte Carlo 
simulation of charm production, on the vertex detector (which was upgraded 
during the run), and on the simulation of the detector stability. We then define 
the split sample variance as the difference between the scaled variance and the 
statistical variance if the former exceeds the latter. The method is described in 
detail in reference [13]. We vary the subresonant states in the Monte Carlo and 
use the variance in the branching ratios as a contribution to the systematic er- 
ror. We investigate the systematic effects based on different fitting procedures 
and we find this contribution to be small. The branching ratio is evaluated 
under various cut selection criteria, and the variance of the results is used as 
an additional systematic error. The systematic effects are then all added in 
quadrature to obtain the final systematic error. Table 1 summarizes the con- 
tributions to the systematic errors for the p^+^^o]j--l+ 7 I+ > l branching ratio. 



The result, 



0.18 ± 0.04 ± 0.05, is summarized in Table 2. 



IV. D+ 



CHANNEL 



This is a challenging channel to reconstruct as we typically only have a de- 
tached silicon track from the production vertex and a K s to indicate a can- 
didate. Several criteria are used to improve the signal over background. Since 
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Table 1 

Summary of the systematic error contributions for Df — ► Kgir~ir + iT + . 



Contribution 


F(D+~*K t >,K-n+n+) 


D Momentum & Run Period 


0.04 


Split Target (in versus out) 


0.03 


Set of Cuts Selection 


0.01 


Fit Variance 


0.01 


Resonant Substructure 


0.01 


Total 


0.05 



any signal was expected to be small the selection criteria are optimized us- 
ing Monte Carlo signal events and sideband background events. The figure of 
merit used was S/y/B and the cuts were chosen sequentially. At each step, 
the S/ \f~B distribution was determined for the full range of each cut. The cut 
which had the highest S/ yf~B was selected and a cut was made more conserva- 
tive than the maximum S/y/B point. The procedure was then repeated until 
no further improvement was possible. 

For the 90% of the K° s decays that occur after the K$ has passed through the 
silicon strip detector, we employ a specialized vertex algorithm to locate the 
Kg7i + vertex. We use the momentum information from the Kg decay and the 
silicon track of the pion to form a candidate Df vector. This vector is inter- 
sected with candidate production vertices which are formed with two other 
silicon tracks. When the D vector is forced to originate at the production ver- 
tex, we can compute a confidence level that the vector formed a vertex 
with the charged daughter. As the type and resolution of Kg is integral to 
finding the Df vertex, the significance of separation, L/<jl, between the pro- 
duction and Df decay vertices were varied according to the K s decay type. 
The L/ol cuts varied from 7-11. This mode also required Iso2 < 2%. 

The normalization channel is the Cabibbo favored Df — > K S K + . The selec- 
tion criteria for this channel (with the exception of particle identification) are 
identical to Df — > Kg7T + . The momentum of the Df and the charged hadron 
in the Df decay must be greater than 45 GeV/c and 12 GeV/c, respectively. 
To reduce the effect of long-lived decays and reinteractions, the proper decay 
time must be less than 2.5 ps with an uncertainty less than 0.12 ps. To help 
separate charm from combinatoric background, a momentum asymmetry cut 



on the two body Df decay was used: 



p(K°)-p{h+) 
p(K»)+p(h+) 



< 0.75. 



For the K + candidate the negative log-likelihood kaon hypothesis, Wk = —2 
ln(kaon likelihood) must be favored over the corresponding pion hypothesis W v 
by W n — Wk > 4 while for the signal mode, the ir + candidate must have Wk — 
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W w > — 1. The first cut serves to dramatically reduce the potentially large 
D + — > Kgiv + background which peaks at the Df mass when reconstructed 
as KgK + while the second cut reduces Df-, — > K$K + background which is 
smaller to begin with and peaks below the Df mass when reconstructed as 

Df - tf°7T+. 

Fitting the Df -> ^tt+ mass plot is complicated by the presence of the large 
L> + — > Kg7r + signal. Since the resolution of the state is relatively poor (a ~ 
13 MeV/c 2 ) there is very little space between the D + and Df peaks to estimate 
the background. The fit used to obtain the central value has five contributions. 
The first contribution is the D + — > K$it + signal which is fit with a distribution 
obtained from smoothing a Monte Carlo sample of reconstructed D + — > Kgir + 
events. The mean and yield are fitted parameters. The second contribution 
is the Df — > Kgir + signal which is also fit with a distribution obtained from 
smoothing a Monte Carlo sample of reconstructed Df — > Kgir + events. In 
this case, the mean is fixed. The third and fourth contributions are reflections 
from Df — + K S K + and D + — > KgK + . The reflection shapes are obtained 
from Monte Carlo samples of generated Dfx — > K®K + events reconstructed 
as Df — > Kgn + . The level is found by taking the same generated events, 
reconstructing them properly, and determining the yield. This Monte Carlo 
yield is then compared to the yield of the data Df — > KgK + and D + — > -Ks-fr' 4 " 
and this factor multiplies the reflection shapes. Finally, the fifth contribution 
is a quadratic polynomial to account for generic combinatorial background. 

The KgK + mass plot is also fit with five contributions. The Df — > i^-ft'" 1 " and 
D + — > -ftlj-fr^ are fit with functions obtained from smoothing reconstructed 
Monte Carlo samples. The masses and yields are fitted in both cases. The 
reflection from D + — > K° s -r + is also obtained from Monte Carlo and fixed 
based on the number of reconstructed D + — > Kg7i + events in data. The fourth 
contribution, a reflection from Df ^K° s K + -n° is allowed in the fit. The shape 
is obtained from Monte Carlo simulation but the level is allowed to vary in 
the fit since the branching ratio is poorly known and we do not have a fully 
reconstructed sample available. As before, the fifth contribution is generic 
combinatoric background which is modeled with a quadratic polynomial. 

From the Kgir + fit shown in Fig. 3 we obtain a Df yield of 113 ± 26 events. 
The K° S K + fit presented in Fig. 3 gives a yield of 777 ± 36 Df events and the 
number of events found for the Df — > KgK + TT° reflection is consistent with 
PDG branching ratios and our efficiency 

The systematic uncertainties are divided into cut variants and fit variants. 
In both cases the systematic uncertainty is obtained from the square root of 
the standard deviation of the values weighted by the individual uncertainty. 
The actual procedure is as follows. For each variant (but not the default), the 
branching ratio BR4 is calculated along with the uncertainty <7j. The average, 
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1.9 1.95 
M(K S K + ) 



1.9 1.95 
M(K s 7t + ) 



Fie. 3. Invariant mass distributions for K° S K + (left) and K%tt + (right). The fits 
are over the entire mass range. Most of the background is modeled by a quadratic 
polynomial. The remaining background is due to reflections and is a different shade. 
The KgK + mode has a large reflection component from Df^>K®K + Tr° below the 
D + peak and a small reflection component from D + — » K® 7r + under the Df peak. 
The 7r + has small reflection contributions below (under) the D + peak from Kgir + 
decays from D + {Df). All signal and reflection shapes come from a Monte Carlo 
simulation. 

weighted by the inverse of the square of the uncertainty, is calculated 



BR 



BRi 



The systematic uncertainty is obtained from the square root of the standard 
deviation which comes from a "weighted" x 2 '- 



a, 



MVS 



\ 



2 BRj — BR 



N-l 



(2) 



where <r is the uncertainty on the default measurement. 



K° S K + sam- 



For each of the cut variants, both the Df — > -R'gvr" 1 " and 
pies are changed the same (with the exception of particle identification cuts). 
The variations are consistent with statistical fluctuations and the systematic 
uncertainty is determined from the standard deviation which is dominated by 
the -D+ — > Kgir + variations. The systematic uncertainty from the cut variant 
is a^ s = 0.010. 



The systematic uncertainty in estimating the yield of — > KgK + events is 
negligible compared to estimating the yield of Df — > Kgir + events. Therefore, 
for the fit variants we vary how the Kg7i + mass plot is fitted. Some of the 
variations include fitting with a Gaussian, allowing the mass and width to 
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float, and fitting only above the D + mass peak. The variation in the Df — > 
Kg7r + yield, again weighted by the uncertainty squared, gives the systematic 
uncertainty. The systematic uncertainty on the yield from the fit variations is 
9.0 events which corresponds to a relative uncertainty of 8.0% and translates 
into a systematic uncertainty on the branching ratio of o{y S = 0.008. Adding 
the cut and fit systematic uncertainties in quadrature gives a total systematic 
uncertainty on the branching ratio of 0.013. 



V. SUMMARY OF RESULTS 



In conclusion we have presented the first evidence of the Cabibbo suppressed 
decay mode Df — > K s 7t~7t + tt + and measured the relative branching ratio of 

r(D+->KU-n+n+) 



-. — - — 0.18 ± 0.04 ± 0.05. A naive expectation for this branching 

1 (VJ —>KgK 7T+7I-+) 1 u 

ratio is tan 2 #c = 0.054. Compared with this expectation the branching ratio 
is more than 3 times larger. One contributing factor is there is more phase 
space available in the D+ — > K s it~it + it + decay than in the D+ — > K s K~tt + tt + 
decay. Another factor is that the Kg in the denominator of the ratio comes 
from a K°. In the numerator the K s may be the result of either a K° or 
a K° decay. Perhaps a better understanding of this ratio would result from 



reporting the ratio 



r(£>+-tf°7r 



Using the branching 



ratio reported in reference [111 for l%i^St + * + j2 = - L : — ± () -° 13 

1 L 1 r(DJ — >KgK 7r+7r+) 

*K%-K-K+-K+) 



we find 



r(r»3 



r(D+^,R'g^-7r+7r+)+r(D+^XgX+7r+7r-) 



0.11. 



We also present evidence for D 
relative to Dt 



■xgi<r+) 



Kg7c^~ and measure its branching fraction 
0.104 ±0.024 ±0.013. This branching 



ratio is also larger than tan 2 8c, but is slightly smaller than predictions [4,5,6] 
which range from 14% to 17%. The results are summarized in Table 2. 

Table 2 

Branching ratios, event yields, and efficiency ratios for modes involving a K s . All 
branching ratios are inclusive of subresonant modes. 



Decay Mode 


Ratio of Events 


Efficiency Ratio 


Branching Ratio 


r(D+^K%ir-ir + iT + ) 
r(D+^K° s K-TT+n+) 
T(D+^K° s n+) 
r(B+-JC°K-+) 


179±36 
763±32 

113±26 
777±36 


1.34 
1.39 


0.18 ±0.04 ±0.05 
0.104 ±0.024 ±0.013 
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